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Abstract 

We present a simple and flexible method of keeping track of the complex phases 
and spin quantization axes for half-spin initial- and final-state Weyl spinors in scat- 
tering amplitudes of Standard Model high energy physics processes. Both cases 
of massless and massive spinors are discussed. The method is demonstrated and 
checked numerically for spin correlations in rf production and decay. Its main ap- 
plication will be in the forthcoming work of combining effects due to multiple photon 
emission (exponentiation) and spin, embodied in the Monte Carlo event generators 
for production and decay of unstable fermions such as the r lepton, t-quark and 
hypothetical new heavy particles. 
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1 Introduction: What is the problem? 

Weyl spinors (Ws) are used in very powerful algorithms to calculate multiple-particle 
and in particular multiple-photon spin amplitudes [0,0. In Ws calculations the complete 
control of the spin-dependent complex phase, i.e. of spin quantization frames, is usually 
neglected. This is not a problem, when one is interested in cross sections for stable 
particles, but it becomes a crucial limitation if one needs to combine them with the spin 
amplitudes for the decay of an almost stable (narrow resonance) fermion. 

The complete spin amplitude for production and decay is given by the formula: 

M= Mf od Mf ca7 , (1) 

A=±l/2 

where A is the projection of the intrinsic spin on the z-axis in a certain fermion rest frame 
(see discussion in the appendix of ref. ||, and also ref. j3J). Let us call this particular 
fermion rest frame the spin quantization reference frame, SQRF. When, in the above 
formula, both production and decay spin amplitudes are obtained from Feynman rules in 
the same calculation, with the same intermediate fermion states u(p, A), the situation is 
simple - we even do not need to know very precisely the SQRFs. The above is not true 
in the many important practical cases such as the production of the r lepton, t-quark 
and hypothetical supersymmetric particles. The reason is the following: the fermion may 
have many decay modes, with complicated multiparticle final states, and complicated 
half-phenomenological decay matrix elements, calculated and implemented separately. 
For instance in the case of r decay, the TAUOLA Monte Carlo || represents by itself 
an independent computational project. In fact every one out of more than twenty decay 
modes of the r represents an independent project both from the computational and from 
the theoretical perspective. 

For these practical reasons, combining production and decay at the density matrix 
level, even if it seems at first more complicated, is worth an effort, because formula ([I]) 
would lead to independent sets of spin amplitudes for distinct combinations of decay 
modes of unstable final-state fermions. The density matrix formalism allows for clean 
modularity in the calculations and in the architecture of the corresponding Monte Carlo 
event generators. 

The price one has to pay (for spin effects not to be destroyed) is that the complete 
combined density matrix for all final-state unstable particles must be provided by the pro- 
duction program with the careful definition of the corresponding SQRFs. That requires, 
in particular, full control of the relative phases of all production spin amplitudes. 

In our paper we shall rely on massless and massive Weyl spinors defined in two papers 
of Kleiss and Stirling (KS) |||7|]. An analogous set of massless spinors is also defined 
in Ref. ||. By supplementing these techniques with rules for complete control of spin 
quantization axes (complex phases) we turn them into a powerful tool for calculations of 
unstable half-spin fermions, useful for a Monte Carlo simulation of production and decay. 
The aim of the present work is to adapt the spin amplitude techniques to the framework 
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established in the works ]|,[||10||> see a ls° 111, 1^1- We shall use it in the forthcoming work 



on exclusive coherent exponentiation |13 



Let us start with the fundamental question of defining properly the spin quantization 
axes (and complex phases) for half-spin Weyl spinors. We present details of the lowest- 
order r-pair production process 

e~Oi,Ai) + e + Go 2 ,A 2 ) -> T'{q u ^ x ) + r + (g 2 ,/i 2 ) (2) 

with spin amplitudes J^x°\ 2 u\u2 wnere A t i>/ i 2 = ±1 denote twice the spin ^-projection of 
the corresponding r lepton in its rest frame, and A, are analogous spin indices of e^. To be 
clear, we will not necessarily always restrict ourselves to Aj and ^ being the conventional 
Jacob- Wick helicities 



The full differential cross section for the combined r-pair production and decay reads 
as follows 

1 1 1 2 \a,b=0 J 

where da prod is the unpolarized differential cross section for the production process only, 
and o?r^ ecay are unpolarized differential partial widths for decays of the t ± . The polarime- 
ter vectors hi are related to the r decay matrix element (see eq. (2.18) of ref. for its 
precise definition)^, and £j are conventional spin polarization vectors of e ± parametrizing 
their spin density matrices p A j = (1 + a ■ e) x \/2. 

All spin effects are embodied in the last term, which includes the full spin density 
matrix related to the production spin amplitudes as follows: 



Rabcd-N 1 ^2 ^A 1 A 2M iM 2 (^A 1 A 2 /Z 1 M 2 )*^A 1 A 1 (7 A2A2 'AlW^ 



N=J2 l- A/i A 1 A2^M2 1 2 , a > b,c,d= 0, 1, 2, 3, 



(4) 



where a k for k = 1, 2, 3 are Pauli matrices and cr° = is just the unit matrix. 

The translation in eq. (fD from spin amplitudes to R a bcd, and the analogous translation 
from decay spin amplitudes to the vectors h c and hj, occurs in the SQRF of r~ for indices 
c, Hx-, [i\ and in the SQRF of r + for indices d, /z 2 , /i 2 . Let us stress again that this 
translation occurs not in some arbitrary r ± rest frame but in precisely the same SQRF 
of t^, where /ij was originally defined as a spin projection on the z-axis! Not only the 
z-axis has to be known, but also the x- and ?/-axesQ. Being in the SQRF makes it also 

1 The vectors hf depend, in general, on momenta of all r decay products. 

2 The change of SQRF by a rotation around the spin quantization z-axis by an angle </> introduces 
the factor exp((«^/2)/ii) (rotation in spinor representation) in the spin amplitudes and consequently an 
i?3(0) rotation in normal space (i.e. adjoint to spinor representation) acting on the index c in R a bcd- 
N.B. a simplest example of Wick rotation is actually i?3(0). 
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legitimate to use Pauli matrices for the translation of spin indices into vector indices in 
eq. 

In the practical calculations that implement the above methodology, see KORALB 
||T0f and KORALZ |16|], a special subprogram TraLor is available, which implements the 
Lorentz transformation from the SQRF[] of each r ± down to the laboratory system. It 
is used to transform the momenta of all decay products into the laboratory frame. 
This arrangement fully guarantees that all spin effects, including spin correlations, are 
properly reproduced in the Monte Carlo event generation. 

Having convinced the reader that the knowledge and the proper use of the SQRF for 
each fermion is very important, let us look into how it is defined. The obvious choice is to 
employ the standard definition of spin states by Jacob and Wick (JW) [14[], the so-called 
helicity states. In ref. || spinors of outgoing r 1 * 1 (and also of e ± beams) were defined using 
a variant of the JW prescription^. In the JW method, states with definite spin projection 
on the z-axis (helicity states) are defined by construction in the SQRF; the full Lorentz 
transformation down to the laboratory frame is used as an inherent construction element 
in the definition of the spin state. In KORALB |),10| the JW spin quantization principle 
was also successfully used to define spin states of e ± and r in the presence of one photon. 

On the other hand, looking at the calculations of refs. ||,|nj it is quite clear that the 
JW method is well suited to lowest order e"e + — > t~t + , but that it is rather complicated 
to deal with already in the case e~e + — > r~r +r ) since the orientation of SQRF systems 
depends on the photon direction even in the soft-photon limit. One definitely needs a 
better method to go beyond the single-photon case. The natural candidates are methods 
based on the Weyl spinors, which have provend to be very successful for the evaluation 
of multiple photon emission spin amplitudes for massless and massive spinors. Generally 
Weyl spinors (Ws) have their own internal definition of the SQRF, substantially different 
from that in the Jacob- Wick method. For example in the KS technique [|5],|7| the spin 
quantization axis in the rest frame of the massive fermion is related to some auxiliary 
light-like vector k involved in the definition of the massive spinor. The typical defini- 
tion of the spinor in the KS method involves a single "brutal" transformation from the 
reference fermion state into the state of definite mass, spin and momentum rather than 
the explicit "soft" rotation/boost method of the JW method. The problem is that in the 
transformation method of the Ws/KS techniques the phase, and therefore the position of 
the x- and y-axes, is lost or undefined. This is not a great problem for most of the practical 
applications of the Ws/KS techniques, such as the calculation of multiple photon emis- 
sion spin amplitudes, which are usually done with the "modulo phase factor" approach 
anyway. This feature of Ws/KS methods inhibits, however, the use of this method for 

3 Pauli matrices play here the role of Clebsch-Gordan coefficients to combine two half-spin representa- 
tions of the rotation group (two indices in the spin density matrix) into spin and 1 angular momentum 
objects, see eq. (A. 10) in ref. pj)j . 

4 The TraLor subprogram actually provides explicit definition of the SQRF frame. 

5 The JW @] recipe was modified in ref. ||, and in fact this choice is closer to that in ref. ||. For 
final fcrmions the y-axis was placed in the reaction plane, while in the JW prescription it is related by 
direct rotation (around the axis perpendicular to the reaction plane) to the laboratory y-axis. 
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unstable fermions. The aim of this paper is to show how to restore the knowledge of the 
complex phase and the SQRF for the KS massive spinors of refs. 0,0]. 

In the following we shall repeat many definitions of refs. in detail. This is made 
necessary by the fact the central question is the control of the spinor phases, and even 
knowledge that this phase is zero (in certain cases) forms an important new result! Our 
final answer for the problem of control of Weyl spinor phases and all three axes of the 
SQRF is very simple - the reader should not be misled by its apparent simplicity, it is an 
important and new supplement of the existing Ws/KS methods. 



2i Basics of the Weyl technique 

For Dirac spinors, wherever an explicit representation is needed, we employ the Weyl 
representation of gamma matrices 
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In the Weyl representation, spinors transform under rotations around the fc-th axis and 
boosts along the j-th follows WA'- 



S(R k ((f>)) = exp 
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S( B j(x)) = exp ( - x 











(6) 



We define eigenstates with definite massless four-momentum p 2 = and chirality A = ±1 

1 



pu\(p) = 0, uj x u x (p) = u\(p), cj a = -(1 + A75), 



(7) 



with the normalization condition u\(p)u\(p) = uo\ p 1 . The four basic massless spinors are 
defined in a certain primary reference frame (PRF), as follows: 



MCt) 









" 

V2 





u-(Ct) 







-V2 



U-(Ci) 






y/2 




(8) 



and they correspond to particles flying along the z-axis: Ct = (1,0,0,1) and Q = 
(1,0,0,-1). It is important to remember that the relative phases in the above set of 
basic vectors are constrained uniquely by the following relations 



u±(C|) = S(i? 2 (7r))u ± (C T ), u + (C t ) = -^u_(Ct), u+(Ci) = jMC), 

whereQ r] = e x = (0, 1, 0, 0), rf = -1, r] ■ ( T = and 77 • C4. = 0. 

6 Our C| corresponds to the vector called ko in ref. M and 77 is denoted there by k\. 



(9) 



4 



3 Massless spinors 



The arbitrary massless spinor of momentum p and chirality A is generated according to 
the KS method , out of the two constant basic spinors of opposite chirality (two out of 
four in eq. (H)) 



u\{p) 



1 



Written explicitly, see also ref. ||, it looks as follows: 

fp 

e i(j> 

«-(p) 



u-a(C), C = Ci- 



(10) 



u+(p) 



p e 










p e 
P 



-i<j> 



p ± =p°± p 3 , p 1 + ip l = p T e l(p . (11) 



The above expressions lead directly to an explicit expression for the inner product of the 
two massless spinors 



s+ipi,Pi) = u + {pi)u-(p 2 ) = \JPiPte t4>1 - \/p s /v 
s_(p!,p 2 ) = «-(pi)w + (p 2 ) = -(s+(pi,p 2 ))*- 



(12) 



It is understood that the above formula is evaluated using four-momenta in the PRF, see 
eq. ([Bp. The same formula can be obtained using the original KS expressionP]: 
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[(pC)(w) - (pv)(q() - ze^C^W 1 



V2K V2gC 

which is applicable in any reference frame, not only in the PRF. 



(13) 



4 Massive spinors 

Spinors for the massive particle with four momentum p (with p 2 = m 2 ) can also be defined 
using a construction similar to that in the massless case: 

u & A ) = ~^=^ (P + m ) U-a(C), v(p,\) = —j=L= (p-m) u A (C). (14) 

As we see, in order to obtain the f-spinor for the antiparticle we have to substitute 
m — > — m and A — > —A, and this is a general rule in the following. Apart from the Dirac 
equation (p — m)u(p, A) = and m)v(p, A) = the above spinors obey the following 
normalization and completeness relations: 

u(p,A)m(p,A) = ^(1 + A 75 *0 (p + m), v(p,A)v(p,A) = -(l + A75tf(jJ-m), (15) 

7 The analogous expression in eq. (3.9) of ref. |(| misses factor 2. For the Levi-Civita tensor we take 
^0123 = 1- 
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where 



1. 



(16) 



is the spin quantization vector, in the fermion rest frame, it points opposite to (, i.e. 
s ~ — (, or in other words the spin quantization axis for every fermion is guided by the 
single and common massless vector (. 

The definition of eq. ([14]) can be rewritten explicitly in terms of massless spinors as 
follows 



m 



u(p, A) = u x (pc) + —, =u_a(C), v(p, A) = w_ A (p c ) 

v 2 K 



m 



V2K 



ua(C), 



(17) 



where 



PC = V - C 



PC 







(18) 



is the light-cone projection of p obtained with the help of the auxiliary vector (. Equa- 
tion (]17D immediately provides us also with the explicit inner product for the massive 
spinors: 



u(px, Xi)u(p 2 , A 2 ) = S(p 1 ,m 1 , Ai,p 2 , m 2 , A 2 ) 
tZ(p x , X 1 )v(p 2 ,X 2 ) = S(px,mi, \!,p 2 , -m 2 , -\ 2 ) 
v(p!, \i)u(p 2 , X 2 ) = S(pi, -mi, -Xi,p 2 , m 2 , A 2 ) 
v(p u X 1 )v(p 2 , A 2 ) = S(p u -mi, -Ai,p 2) -w 2 , -A 2 ) 



(19) 



where 



5(pi,mi, Ai,p 2 ,m 2 , A 2 ) = 5 Al _ A2 s Al (pi c ,p 2c ) + 5 Al)A2 I "11^/^77^ + m 2\lT^- ) • (20) 



2CPi 



2CP2 



For the sake of the subsequent discussion, let us also write the explicit basic massive 
spinors in the particular case when the fermion rest frame coincides with the PRF (in 
which C = (1,0, 0,-1)): 



u{p, +) 



u{p, 



v(p, +) 





-1 

1 



v(p, -) 



1 


-1 




(21) 



We have omitted the trivial normalization factor ym. The convention for relative phases 
is inherited from the underlying massless spinors @. In particular the important phase 
relations 



S(R 2 {-K))u(p,+) = u(p, 



S(R 2 {7r))v(p,+) = v(p,-). 



(22) 
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For ultrarelativistic p° — > oo, or for an almost massless fermion m — > 0, the differ- 
ence between p and p^ is negligible and the massive spinor becomes identical with the 
corresponding massless spinor: 

lim u(p, A) = lim v (p, — A) = U\(pA, 

m— »0 m^O (23) 

lim S(pi,m 1 , \i,p2,m 2 , A 2 ) = <5 Al -a 2 sa 1 (pi c ,P2 C ) 

The validity of the above limits is restricted by the condition p ■ ( ^> m, i.e. they fail 
in the special case of p parallel to ( and such cases have to be treated one by one, with 
special care. In such a case, in eq. (JlTf) , the second term, proportional to mass, survives 
and dominates. The same phenomenon occurs in the inner product of eq. (|2"0D . Thus, we 
have to remember that the naive procedure of taking the massless fermion limit by means 
of straightforward omission of all terms proportional to mass can be dangerous and the 
limiting procedure has to be done very carefully^. 

5 Again the central question and the GPS answer 

We now come to the central question of interpreting the massive fermion spinors of 
eq. (0), the spin vector of eq. (|T6|), and finding out what the SQRF for these spinors 
is. The states u(p, ±) and v(p, ±) represent pure quantum mechanical states, which we 
intend to use in the Feynman rules to calculate scattering matrix elements. As we see 
for instance from eq. (|16D , u(p, +) and u(p, — ) have spin projections +1/2 and —1/2 on 
the vector s in the fermion rest frame. Thus the vector s, which is guided by (, deter- 
mines the z-axis of the SQRF. This is the easy part of determining the SQRF. The more 
difficult question is: Where are the x- and y-axes of the SQRF corresponding to u(p, +) 
and u(p,—) defined in eq. (|14D ? Is this question meaningful at all? Yes, it is meaningful 
because a rotation Rz{(j)) around the z-axis introduces a A-dependent phase in u'(p, A). In 
fact one possible method of finding the x- and y-a.xes of the SQRF is to inspect u(p, ±) in 
the rest frame of the fermion and to adjust such that the relative phase of u'(p, +) and 
u'(p, — ) is zero, i.e. such that the relation S(R2(7i))u'(p, +) = u'(p, — ) holds. This rule is, 
however, not very practical - we would like to have at our disposal more practical, easy 
to apply, rules, which determine the position of the same x- and y-axes of the SQRF for 



both (A = ±) spinors defined in eq. (|14l) . 



The rules for determining all three spin quantization axes for u(p, ±) and v (p, ±) of 



eq. (|14]) are the following: 

• In the rest frame of the fermion, take the z-axis along 

• Place the x-axis in the plane defined by the z-axis from the previous point and the 
vector rf, in the same half-plane as ff. 

• With the y-axis, complete the right-handed system of coordinates. 



8 On the other hand, keeping mass terms does not necessarily mean correct results in numerical cal- 
culations, because of machine rounding errors. 
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We call the above rules the Global Positioning of Spin scheme or in short the GPS rules, 
and we shall call the GPS frame the SQRF determined by the GPS rules. We give here 
the formal proofs of the above rules. 

Proof: Since all massive spinors are defined in terms of two basic massless spinors 
U +(C|) an d U_(Cj,), it is sufficient to consider underlying massless spinors and to prove 
that after transforming them from the PRF to the GPS frame they look the same as in 
eq. (H), up to a common phase factor. In other words the two basic massless spinors 
u + (Cj.) and u_(Cj,) get regenerated in the GPS frame. Without actually doing explicitly 
the transformation from the PRF to the GPS frame we know that because of the Dirac 
equation, chirality deefinition (see eq. (|7|) ), and because ( = Q = const (1,0,0, —1) the 
resulting spinors have to look as follows 



U-(C|) 



(24) 



where C2 and C3 are unknown complex constants. We shall now exploit two facts: (a) in 
the local GPS frame, ff lies in z-x plane (with positive x-component) and (b) the four- 
vector i] obeys 77 • (4 = and rf = — 1. The above implies that 77 = (—a, 1, 0, a), where a 
is a real constant, and therefore 



jf 



The above completes the proof because we may exploit the relation 
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-2a 


-2a 


-1 








-1 












(25) 



u+(C|)=tfu_(C|) 



(26) 



see eq. (^) , finding out that 03 = 02, and therefore we do really get back the same looking 
basic spinors in the GPS frame as in the original PRF, up to a spin-independent phase 
factor c = c 3 = C2 = rexp(z0). In fact one may advance even one step further and prove 
that the actual overall phase <ft is exactly zero. We sketch a formal proof of this conjecture 
in the following. 

Proof: Let us observe that the spinor transformations S(Bi), S(B 3 ) and 5 , (i?2), i-e. 
transformations induced by boosts in the x-z plane, and the rotation around the y-axis 
are purely real, while the S(R 3 ) induced by the rotation around the z-axis is complex and 
diagonal, see eqs. @. We shall now analyse the complete string of the Lorentz trans- 
formations connecting the PRF and the GPS frames of a given fermion with momentum 
p. The first rotation i?3(</> p ) puts p in the x-z plane and introduces the following phases: 
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Ua(C|) — > exp(i0 p /2A)UA(Cj.)- Next, the rotation R 2 and boost B 3 and again the rotation 
i?2 transform us to the frame where p = (m, 0,0,0) and Q = (1,0,0,-1). Since the 
corresponding transformations of spinors are real, we know that the resulting final Ua(Cj.) 
preserves its complex phase! We are not yet in the GPS frame - the last transformation 
needed is the rotation -R3 (</>') which brings if into the x-z plane. The resulting net phase 
is exp((i/2)(0 p + 0')A). Since we already know from the former proof of the GPS rules 
that there is no A dependence in the total phase, the only possible solution is that the 
total change of the spinor phases is zero, <p p + <p' = 0, modulo 4wr for some integer n, so 
that exp((i/2)(0 p + (p')X) = 1 as desired. (It is also possible to give the argument, based 
on geometry, that P = —0'.) We have also checked numerically, using spinor transforma- 
tions of eq. ([]), that the series of the above Lorentz transformation really regenerates the 
basic spinors u±(C) in the GPS frame. For massive spinors, we have checked numerically 
that the massive spinors defined in eq. ( [I?D in the laboratory frame, when transformed 
properly to the GPS frame, become identical to those of eq. fl21~|) ! 

Summarizing: the important property of the GPS rules is that the basic spinors defined 
in eq. fll4|), when transformed to the GPS rest frame of the fermion, obtain the same form 
(are regenerated), up to a real spin-independent, normalization constant, as the original 
ones, defined in eq. fl2~T|). The Lorentz transformation from any Lorentz frame, for instance 
from the laboratory reference system, to the SQRF/GPS frame of every massive fermion 
is therefore uniquely defined. Thanks to the GPS rules, we are finally in a situation as 
comfortable as in the Jacob- Wick method. The difference is however clear: in the JW 
method, the transformation from the laboratory frame to the SQRF is the essence and 
the starting point of the method - it is tied up to the fermion four-momentum in the 
laboratory system, and the axes of that frame. In the GPS method the SQRF/GPS 
frame and the Lorentz transformation from the laboratory to the SQRF/GPS frame are 
deduced with the help of the "backward engineering" procedure guided by the auxiliary 
vectors ( and 77 (common to all fermions), without any direct reference to the laboratory 
system. 

6 GPS at work, r-pair spin correlations 

In the following we shall show how to apply the GPS rules in order to find out the Wick 
rotation || responsible for change of spin quantization axes between the JW and KS 
conventions for fermion spinors, for the spin amplitudes of the e~e + — > t~t + process. 

Let us use the variant of JW helicity states defined in ref. and refer to it as JW2. 
The JW2 spin states are defined (quantized) in the r~ rest frame with axes: ~z = — t+ 
and y = r+ x e~/|r + x e~|, and in the r + rest frame with axes: z = —t~ and y = 
—t~ x e + /|r~ x e + |. We denote as / the unit three-vector parallel to the momentum of 
a fermion /. The other (third) axis is always picked in such a way that the right-handed 
reference frame (x, y, z) is formed. With the above choice of quantization axes the double 
spin density matrix of the r pair (see formula (2.6) of ref. 0, used also in the KORALB 



9 



+++++++++++++++++++++++++++++++++++++++ Born spin amplitudes +++++++++++++++++++++++++++++++ 



£111 
£ 1-1 1 
£-1 1 1 
£-1-1 1 
.000000 
.000000 
.000000 
.000000 



o 
15 
1} 
15 



£ 1 1-1 
£ 1-1-1 
£-1 1-1 
£-1-1-1 



o 
o 

15 

15 



£11 
£ 1-1 
£-1 1 
£-1-1 



1-15 
1-15 
1-15 
1-15 



£ 1 1-1-15 
£ 1-1-1-15 
£-1 1-1-15 
£-1-1-1-15 



.000000] [ 

-.454758] [ 

.454758] [ 

.000000] [ 



.000000 
-.696227 
-.280666 

.000000 



.000000] 
-.001005] 
.001005] 
.000000] 



[ .000000 

[ -.280666 

[ -.696227 

[ .000000 



.000000] [ 

-.001005] [ 

.001005] [ 

.000000] [ 



.000000 .000000] 

.000000 .454758] 

.000000 -.454758] 

.000000 .000000] 



0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 

##################################### R-matrix GPS fram ####tt###tt####################tt###tt### 

[ 1.000000 .000000] [ .000000 .000000] [ .000000 .000000] [ .000000 .000000] 
[ .000000 .000000] [ -.023328 .000000] [ .000855 .000000] C .001871 .000000] 
[ .000000 .000000] [ -.000855 .000000] [ .823264 .000000] [ -.386811 .000000] 
[ .000000 .000000] [ .001871 .000000] [ .386811 .000000] [ -.153408 .000000] 

##########tt###tt#tt###tt###tt##########tt# R-matrix JacHick ####tt###tt#tt###tt###tt##########tt###tt#tt# 

n 1.000000 .oooooo] : .000000 .000000] : .000000 .000000] : .000000 .000000] 

[ .000000 .000000] [ .199012 .000000] [ .000000 .000000] [ -.375031 .000000] 
[ .000000 .000000] [ .000000 .000000] [ .023418 .000000] [ .000000 .000000] 
[ .000000 .000000] [ .375031 .000000] [ .000000 .000000] [ -.824406 .000000] 



******************* R analytical 



******************* 



C 1.000000000000] [ 

[ .000000000000] c 

[ .000000000000] [ 

[ .000000000000] [ 



.000000000000] [ 

.199012228374] [ 

.000000000000] [ 

.375031393988] [ 



.000000000000] [ 

.000000000000] [ 

.023418353587] [ 

.000000000000] [ 



******************* R Difference 



******************* 



[ .000000000000] [ 
[ .000000000000] [ 
[ .000000000000] [ 
[ .000000000000] [ 
cos (theta ) ,sin (theta )= 



.000000000000] [ 

.000000000000] [ 

.000000000000] [ 

.000000000000] [ 
-.884680271149 



.000000000000] 
.000000000000] 
.000000000000] 
.000000000000] 
.466198260229 



amfin,sqrt(s)= 



1.777000000000 4.000000000000 



.000000000000] 
.375031393988] 
.000000000000] 
.824406125213] 



.000000000000] 
.000000000000] 
.000000000000] 
.000000000000] 



Figure 1: Numerical illustration of Wick rotation for one r-pair production event at y^s =4 GeV, 
m T = 1.777 GeV, cos# = —0.884680271149. In the upper part we print the complex spin amplitudes 
■M-^ , with the numbering style {AiA2/^i^2} shown explicitly. Then, we show the Roocd-, c,d = 
0...4, correlation matrix (imaginary part equal zero) as calculated in the GPS frame, using eqs. (Q) 
and (|28|) , and after a (double) Wick rotation to JW2 r rest frames. On the other hand, the result 
marked "analytical" is obtained from eq. (p7|), also in the JW2 frames. The difference of the two 
results in the JW2 frame is shown to be precisely zero. 



Monte Carlo [10]) takes the form: 



JW2\ 
OOcd J 



J\f 



1 + c 2 + M 2 s 2 








[l + M 2 )t 


-2Mcs 



(1 






M 2 ) 






2Mcs 




-1 - c 2 + M 2 s 2 



(27) 



where c = cos#, s = sin 9, 9 is the scattering angle between e~ and r~ , M = 2m T /s 1 ^ 2 
and M — 1/(1 + c 2 + M 2 s 2 ). In order to simplify the discussion we include only s-channel 
photon exchange (the parity- violating Z exchange is neglected). We also neglected the 
electron mass. To be more precise, in ref. M the spin quantization axes of r are slightly 
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different from JW2; let us call them the KB choice. The KB spin states are defined 
(quantized) in the r~ rest frame with axes z = t+ and x = r + x e~/\r + x e~\ and in 
the r + rest frame with axes z = — t~ and x — t~ x e+/|r _ x e+|. The transition from 
KB to JW2 costs additional trivial Wick rotation R s (n/2) for r + and R 2 { , k)Rz{'^ /^) for 
r~. (These rotations are merely permutations and reflection of axes.) This is what was 
actually done in order to get (analytically) our eq. (|27l) from eq. (2.6) of ref. §|. 

On the other hand, we calculate spin amplitudes for the same process with the basic 
KS spinors of eq. fllTD . Using the Chisholm identity (treating electron spinors as massless) 
and the inner product of eq. (PD|), we obtain the following result: 

MV Xix£i^ = C ^i.-Aa [ U(qi, Hi)u{p!, Ai) V(p 2 , \ 2 )v(q 2 , /i 2 ) 

+u(q u fii)v(p2, -A 2 ) u(pi, -Xi)v(q 2 , /i 2 )] ^ 
= C 5 Xl -\ 2 [ S(q 1 ,m T ,jj 1 ,p 1 ,0,\ 1 ) S(p 2 ,0, -\ 2 ,q 2 , -m T , -// 2 ) 
+S(q 1 ,m T ,ii 1 ,p 2 ,0, A 2 ) S(pi,0, -Xi,q 2 , -m T , -fJ, 2 )], 

where C is an overall normalization constant, unimportant for our discussion. As com- 
pared to the previous result, we lack insight into the analytical structure of the above 
formula; it is simply meant to be evaluated numerically. We are also unable to show ana- 
lytical expressions for -R^f - the translation from M^^^ ^° R-abcd using ec L- (ED has to 
be done numerically. In view of the nice simplicity of the JW/KB helicity amplitudes (see 
eg. eq. (2.3) of ref. |J), one may wonder why bother about all the complications of the 
KS/GPS scheme. The answer is well known: the KS spinors (now with GPS supplement) 
will be very powerful and useful once we consider multiple photon emission from fermions. 
In this case, the Jacob- Wick helicities are no longer so simple and useful, as they were 
really invented for the Born amplitudes. 

The following exercise would demonstrate that our GPS rules really work: If we cal- 
culate Rq£.j numerically and if we really know precisely (from the GPS rules) in which 

reference frames R^Jl is defined, then we should be able by the appropriate two non- 
trivial Wick rotations to get exactly the very simple R J ™J of eq. fl27|). The two Wick 
rotations are in this case just ordinary rotations in two r rest frames, corresponding to 
transformations from the GPS to the JW2 rest frames. In the following, see fig. 1, we 
show the numerical results, demonstrating precisely the above operation done numerically. 
Indeed, we see that due to the appropriate Wick rotations the i?-matrices become identi- 
cal. This demonstrates clearly that in fact our GPS supplement to the KS spinor method 
works in the practical application. The above exercise is just a small warm-up before the 
forthcoming paper in which we apply the GPS rules to a full-scale calculation with 
multiple photon emissions (with exponentiation) and with decaying fermions. 

T Conclusions 

In the present paper we have established a method to control the complex phase and 
spin quantization reference frame for the Weyl spinor technique (Kleiss-Stirling type) for 
massive and massless spinors. The method applies to multiphoton emission. The present 
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work opens the way to extending the work of refs. |5|,|,[T0[ (where the spin effects, photon 
emission and Monte Carlo event generation were successfully combined for the first time) 
to the case of the multiple photon emission. This will be elaborated in a forthcoming 
publication JT3 . 
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